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ABSTRACT 

We track the histories of massive clusters of galaxies fornied within a cosmological 
hydrodynamic simulation. Specifically, we track the time evolution of the energy in 
random bulk motions of the intracluster medium and X-ray measures of cluster struc- 
ture and their relationship to cluster mergers. We aim to assess the viability of the 
turbulent re-acceleration model for the generation of giant radio halos by comparing 
the level of turbulent kinetic energy in simulated clusters with the observed properties 
of radio halo clusters, giving particular attention to the association of radio halos to 
clusters with disturbed X-ray structures. The evolution of X-ray cluster structure and 
turbulence kinetic energy, k, in simulations can then inform us about the expected 
lifetime of radio halos and the fraction of clusters as a function of redshift expected 
to host them. We find strong statistical correlation of disturbed structure measures 
and the presence of enhancements in k. Specifically, quantitatively "disturbed" , radio 
halo-like X-ray morphology in our sample indicates a 92% chance of the cluster in 
question having k elevated to more than twice its minimum value over the cluster's 
life. The typical lifetime of episodes of elevated turbulence is on the order of 1 Gyr, 
though these periods can last 5 Gyrs or more. This variation reflects the wide range of 
cluster histories; while some clusters undergo complex and repeated mergers spending 
a majority of their time in elevated k states, other clusters are relaxed over nearly 
their entire history. We do not find a bimodal relationship between cluster X-ray lu- 
minosity and the total energy in turbulence that might account directly for a bimodal 
Lx — PiAGHz relation. However, our result may be consistent with the observed bi- 
modality, as here we are not including a full treatment of cosmic rays sources and 
magnetic fields. 

Key words: cosmology: theory-galaxies:clusters:general-cosmology:observations- 
hydrodynamics-methods:numerical- X-rays: galaxies: clusters-galaxies: clusters: intr- 
acluster medium 



1 INTRODUCTION 



Galaxy clusters have long been known as radio synchrotron 
sources. Given that observations show that the intraclus- 
ter medium is typical l y magnetized at the level (see 

ICarilli fc Taylor! |2002| : iGovoni fc Ferettilliooi i. what is re- 
quired for the radio emission is a nonthermal population of 
accelerated electrons. Arguments for various scenarios re- 
sponsible for cluster radio emission therefore hinge on the 
method of particle acceleration, using the observed radio 
and X-ray properties of galaxy clusters as discriminatory 
evidence. Radio properties of the sources in galaxy clusters 
vary, though in general the types of sources can be separated 
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into two basic categories, sources associated with individual 
galaxies and those associated with the diffuse intracluster 
medium (ICM). Those associated to individual galaxies in- 
clude the jet and lob e sources connected to active galac- 
tic nu clei (AGN) (e.g.. iHardcastle et al.|[2003l : iNulsen et all 
The diffuse sources can be further divided into two 
broad classes, radio relics and radio halos. Radio relics in- 
clude both steep spectrum re-energized AGN jet lobes, as 
well as polarized sources likely associated with recently ac- 
celerated cosmic-ray (CR) electrons at merger shocks (e.g., 
IClarke fc EnsslinI l2006l : Ivan Weeren et al.l I2OI0I ') . These ex- 
tended, elongated features are steep spectrum sources with 
curved shapes, spatially distinct and strongly argued as co- 
incident with sh ocks from both observations and numeri- 
cal simulations l|Roettiger et al. 1 1 19991 : H assim et al.l I2OOII : 
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iBagchi et al.ll2006l : [SkrUman et al.ll201Cl ). Radio halos (RHs), 
by contrast, are Mpc-scale, diffuse, unpoiarized, steep spec- 
trum radio emitting regions typicaily centered on tiie X- 
ray center of tlie IC M, witli a spatial distribution similar 
to the X-ray e mission (Liang et al. 2000:' Feretti et al.ll20oil : 
iGiacintucciet al. 2009; Macario ct al. 201^. 

The origin of RHs is less certain than the ra- 
dio re l ics, t hough it was recognized a decade ago (e.g. 
iBuotd I2001I ). and observational evidence now continues 
to s uggest that they are associated with merging clus- 
ters (ISarazinll2004l: iGovoni et al.ll20o3 : iBrunetti et"al]|2009l : 



Cassano et al.l 



Cassano et al.l 



For example, the recent work of 
uses radio observations (with GMRT) 
and Chandra observations of an X-ray selected sample of 
galaxy clusters to make a strong case for a division in the 
observed X-ray morphological properties of RH and non- 
RH galaxy clusters. In brief, clusters that are more dis- 
turbed in the X-ray power ratios, centroid shift measures 
and have lower concentration parameters are those that tend 
to host radio halos. Those with relaxed structure measures 
and higher concentrations tend to not have detectable large- 
scale radio emission. 

Recent observational and simulation work argues that 
RHs are created through merger-induced turbulent re- 
acceleration of an existing population of CRs in the ICM. 
The argument for turbulent re-acceleration of a pre-existing 
CR electron population as the cause of RHs appeals to 
a combination of observational properties. First, the radio 
emission is broadly distributed over a relatively large vol- 
ume, roughly centered on the X-ray centroid of the clus- 
ter, and more or less spatially uniform (for a recent review, 
see [Ferrari ct al. 2008). Given the short synchrotron life- 
times of the CR electrons responsible for the emission, they 
must be accelerated m situ, ruling out acceleration from 
large individual shoc k structures on morphological grounds 
(jSrunetti et al.l [20011 ). Turbulence provides a source of en- 
ergy that is distributed relatively evenly over the cluster vol- 
ume. Second, this model naturally explains the connection 
of the presence of RHs to cluster mergers . Though secon d 
order Fermi processes are usually invoked l|Petrosianll200lf ). 
there are legitimate questions about the exact method by 
which the turbulent energy can be converted into acceler- 
ated particles. In addition, a key prediction of the turbulent 
re-acceleration model is th e expected pre sence of very steep 
spectrum RHs (q > 1.6) jCassan 0I I2OO9I '). Two recent RHs 
have been fo und to host so c alled ultra steep spectrum radio 
halos , A521 jBrunetti et al.l i2008') and A697 (jMacario et all 
I2OIOI ). These discoveries, though small in number so far, fit 
the turbulent re-acceleration model. 

Another explanation proposed for the origin of RHs is 
hadronic production. This model invokes the production of 
secondary relativistic electrons generated by the collisions of 
relativistic protons and the thermal population of atomic nu- 
clei in the ICM. In the hadronic scenario, the secondary elec- 
trons are responsible f or the observed radio emission (e.| 
iPfrommer etai]|2008l : [keshet fc Loebl I2OI0I : lKeshetll201C 
and the proponents argue that the hadronic model more 
naturally explains the correlation of the radio emission with 
the X-ray emitting ICM (i.e. the target population for p-p 
collisions) and the presence of both RHs and radio relics 
in clusters. Merger induced enhancement of cluster mag- 
netic fields is invoked to explain the RH-merger connection 



iKushnir et~alll2009l : Ikeshet fc Loebll2010l 'l. However, recent 
7-ray upper limits from Fermi observations of galaxy clus- 
ters are beginning to constrain the capability of hadroni c 
models to account for RHs fsee lJeltema &: Profumci|l201ll '). 
as the hadronic processes should also generate 7-rays from 
neutral pion decay. The recent discovery of ultra steep spec- 
trum radio sources described above is also not predicted 
in the hadronic scenario given the energetic considerations 
required by the 7-ray upper limits. These and other con- 
straints suggested by observations of RHs, such as the spec- 
tral a nd spatial properties of the radio emission (see lBrunettil 
|2004 ). argue against the hadronic model. 

One more piece of observational evidence in the case 
of radio halos is the relationship between the bulk X-ray 
and radio properties of massive galaxy clusters. The radio 
luminosity at l.AGHz and the X-ray luminosity of massive 
radio halo galaxy clusters exhibit a strong correlation for 
many objects. However, there apparently exists a bimodal 
distribution in the L^c — Piaghz plane, where many clus- 
ters have no radio detections, only upper limits in Piaghz 
that lie well below the scaling relation. These objects are of 
similar mass and X-ray luminosity to those that clearly host 
radio halos and are on the scaling relation. There is a signifi- 
cant gap between the measured radio luminosity of the RHs 
and the upper limits of the non-RH objects, which suggests 
that the transition between the two p hases (radio halo ver - 
sus non-RH) must be relatively quick (jBrunetti et al.|[2009l ). 
While there were initial concerns that the bimodality of the 
Lx — PiAGHz relation might be an ar tifact of selection eflfects 
fe.g.. iRudnick fc Lemmermanll2009h . the depth of the radio 
data for the current samples appears to address that issue. 
The presence of a strong correlation of radio power with 
X-ray luminosity for RH clusters, the bimodality of the re- 
lation, and the inferred lifetime and transition times for RHs 
have been used t o argue for or again s t models of radio hal o 
production fe.g. iKushnir et al.|[2009l : IBrunetti et al.1 120091 '). 
In this case, numerical simulations provide a window into 
the evolution of the relevant physical properties of the ICM. 

In this paper, we study the evolution over time of mas- 
sive clusters formed within a cosmological volume simula- 
tion, which allows us to track cluster mergers, the evolution 
of cluster structure, and the presence and evolution of tur- 
bulence. Using simulations, we can investigate the frequency 
and timescale over which clusters exhibit elevated levels of 
turbulence and X-ray structures similar to those observed 
for RH clusters. The use of cosmological simulations turns 
out to be critical to accurately capture cluster merger his- 
tories. In particular, we attempt to understand the connec- 
tion between the X-ray morphology of galaxy clusters and 
the putative energy source for CR acceleration in the ICM: 
turbulence. We do not attempt to determine the method 
of particle acceleration, we merely discuss the probability of 
the turbulent re-acceleration model of RHs given the correla- 
tions between the physical properties of the simulated clus- 
ters and their observed morphological characteristics. We 
study both the "snapshot" correlations for our clusters be- 
tween the X-ray structure measures (power ratios and cen- 
troid shifts) and the turbulence kinetic energy and the time 
evolution of these properties. 

Section 2 of this paper describes the methodology of the 
analysis we undertake, including the details of the numerical 
simulations and analysis and the method of calculating the 
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X-ray structure measures from the synthetic observations. 
Section 3 describes the distribution of simulated galaxy clus- 
ter structures compared to observations, the time evolution 
of both the structure measures and turbulence kinetic en- 
ergy, and compares how the turbulence kinetic energy cor- 
relates with the observed properties of radio halo clusters. 
Section 4 gives estimates for lifetimes and duty cycles for 
radio halos based on the typical duration of episodes of ele- 
vated turbulence kinetic energy and the fraction of clusters 
in these states. Section 5 summarizes and discusses the im- 
plications of our results. 



2 METHODOLOGY 
2.1 Simulations 

For this study we use a numerical cosmological 
hydrodynamic/N-body simulation generated using the 
Enzc0 code, a publicly available adaptive mesh refinement 
(AMR) cos mology code developed by Greg Bryan and 
colleagues jBrvan fc Norman ' l997al lbl: iNorman fc Brvaal 
Il999l : lO'Shea et al.. ,2004 , ,2005 ). The specifics of the Enzo 
code are described in detail in these papers (and references 
therein) . 

This simulation is set up as follows. We initialize our 
calculation at z = 99 assuming a cosmological model with 
Qm = 0.27, Qb = 0.04, Ha = 0.73, h = 0.7 (in units of 
100 k m/s/Mpc), as = 0.9, and using an lEisenstein Sz Hul 
l|l999 !) power spectrum with a spectral index of ris = 0.97. 
The simulation is of a volume of the universe 128 h"^ Mpc 
(comoving) on a side with a 256^ root grid. The dark matter 
particle mass is 9.1 X 10^ h"^ Mq. The simulation was then 
evolved to z = with a maximum of 5 levels of adaptive 
mesh refinement (a maximum spatial resolution of 15.6 h^^ 
comoving kpc) , refining on dark matter and baryon overden- 
sities of 8.0. 

This simulation includes a prescription for radiative 
cooling of the gas using non-equilibrium cooling and chem- 
istry for H and He, and a spatially uniform but time- 
varying Haardt-Madau ultraviolet radiation background. In 
this case, we have not included metal line cooling as has 
been done in previous calculations, since in this study we 
are interested in the morphological properties of the clus- 
ters. Many investigators have documented the over-cooling 
problems associated with including metal cooling in cosmo- 
logical simulations, and in previous work we have shown 
that this over-cooling leads to a severe mismatch of struc- 
ture measures between the observed and simulated cluster 
images. 

We perform all our post-proce ssing analysis for this 
work using the YT0 analysis toolkit l|Turk et al.ll201ll '). For 
each of 132 simulation outputs equally spaced in time {St = 
0.055 Gyr) in the redshift i nterval ^ z ^0.9 , we ran the 
HOP halo-finding algorithm (|Eisenstein fc Hutlil998l ') on the 
dark matter particle distribution to produce a dark matter 
halo catalog. For each halo catalog in each time interval, we 
create spherically averaged radial profiles of a set of physi- 
cal properties (including density and mass) to calculate the 

^ http://cnzo.googlecode.com/ 
^ http://yt.cnzotools.org/ 



value of r2oo and M200 for the identified halos. M200 refers 
to the total mass inside a radius of r2oo , the radius at which 
the overdensity average inside the sphere centered on the 
cluster is 200 times the critical density. We then took all the 
halos at z=0 with M200 3 x lO^^M© (a relatively arbitrary 
cutoff) as our sample. We choose the high mass objects so 
that we can ma ke a reasonable comp arison to the very high 
mass objects of lCassano et al.l (|2010l ). and to limit our sam- 
ple to the best resolved objects in the simulation. This cut 
results in a sample at z=0 of 16 galaxy clusters. 

Following the identification of the sample, we gener- 
ated a set of synthetic 0.3-8.0 keV X-ray i mages of each 
clust er generated from the Cloudj0 software (|Ferland et al.l 
|l998i). We then make a corresponding image of each of these 
clusters from each of the 132 time outputs described above. 
Each initial image represents a projected area of 8 by 8 Mpc 
around each cluster in order to study both the cluster and 
the surrounding structures. The X-ray images are output in 
FITS format with the correct angular scale for their redshift, 
and with the flux modified by the distance. Therefore, we 
are left with a series of synthetic X-ray surface brightness 
images that we can analyze in a similar way to observed im- 
ages of galaxy clusters. Each series represents a time history 
of each of our 16 clusters. Additionally, we have calculated 
the bulk properties of each cluster and its most massive pre- 
decessors throughout this time history. 

2.2 Structure Measures and Turbulence 

2.2.1 Structure Measures 

The association of radio halos with clusters undergoing 
mergers, known visually for quite some time, has been es- 
tablished quantitatively through the presence of significant 
substructure or disturbed cluster structure in the X-ray sur- 
face brightness distribution (e.g. 'Buote"200lVCassan o et al.l 
2010). In particular, Cassano et al. (2010) use three struc- 
ture measures based on X-ray images of the clusters in their 
sample: the centroid shift w, t he third order power ratio 
P3/P0, and the concentration c. ICassano et al.l (|201(]| ) show 
that giant radio halo clusters can be efficiently separated 
from non-radio halo clusters and radio mini-halos in the 
plane defined by any two of these structure measures (i.e. w 
vs. P3/P0, c vs. w, or c vs. Ps/Po). 

Here we will employ the former two statistics, the cen- 
troid shift and P3/P0, the calculations of which have been 
discussed in detail in prev ious papers (|jeltema et al.l [20051 . 
l2008l : iHallman et al.|[2O10l '). In brief, the centroid shift mea- 
sures the variation (here standard deviation) of the position 
of the centroid of the X-ray surface brightness of a cluster 
relative to the X-ray peak in apertures of increasing radius, 
and w is normalized relative to the largest aperture consid- 
ered. The power ratios are based on the multipole moments 
of the X-ray surface brightness in an aperture of a given size; 
specifically they are proportional to the sum o f the squares of 
the m oments divided by the total cluster flux (|Buote fc Tsail 
I1995II . P3/P0 is the third order power ratio and is sensitive 
to deviations f rom m irror symmetry. For consistency with 
ICassano et al.l (|201(]| ). we employ the same overall aperture 
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size of 500 kpc for both w and P3/P0 as well as the same 
aperture spacing of 25 kpc for the centroid shift (see for ex- 
ample equations 1-5 in lCassano et al.l (j2010l )l. In the case of 
the centroid shift, we remove the central 50 kpc surround- 
ing the X-ray peak from consideration. The centroid shift 
w and the power ratios within an aperture of 500 kpc are 
calculated for every time output over the history of our 16 
simulated clusters. 

We do not use the c o ncent ration parameter, which is 
defined by ICassano et all (|2010l ) as the ratio of the X-ray 
surface brightness within the central 100 kpc to the surface 
brightness within 500 kpc. With this definition, the typical 
cluster concentration evolves significantly with redshift as 
a fixed 100 kpc aperture encloses more or less of the clus- 
ter volume. This evolution makes a fixed cut for relaxed or 
disturbed clusters impossible. A different definition of con- 
centration based on overdensity radii would provide a bet- 
ter criterion for clusters spanning a large redshift range but 
would not all ow a direct comparison to the observed radio 
halo sample of lCassano et al. We note, however, that 

in the redshift range spanned by the ICassano" et al.1 (|2010l ) 
sample, we find a similar median concentration for our simu- 
lated clusters to that of the observed clusters once the X-ray 
peak is removed (as was done for the centroid shifts). 

2.2.2 Turbulence Kinetic Energy 

To evaluate the connection between the morphological struc- 
ture measures and the physical state of the cluster, we have 
made a calculation of the mass-averaged turbulence kinetic 
energy (TKE) for each cluster inside the same radius (500 
kpc) that we have calculated the structure measures. Here 
we use a commo n definition of the turbulence kinetic energy 
density (see e.g., IChoudhuri|[l998l ') 

fc = 0.5(^ + Mf + «!), (1) 

where Ux,Uy, are the orthogonal components of the pe- 
culiar velocity of the gas in each grid zone (subtracting the 
bulk halo velocity), and k is the specific kinetic energy. The 
bulk halo velocity is calculated as the mass-weighted aver- 
age velocity of all dark matter particles identified by the halo 
finding algorithm as halo members. For our analysis we cal- 
culate a mass-averaged value of fc in a sphere with r=500 kpc 
around each cluster center. We note that we are not truly 
calculating the energy of the full turbulent cascade. We can 
calculate a physical Reynolds number for flows of character- 
istic length scale equal to the minimum grid scale (for grid 
elements at the peak spatial resolution) of the simulation. 
Our grid size at peak spatial resolution and the character- 
istic parameters of the ICM flows can be combined into a 
dimensionless Reynolds number 

Re = vL/u, (2) 

where v is the fluid velocity, L is a characteristic length, and 
V is the kinematic viscosity of the fluid. We can estimate 
the values for v and L to be roughly the sound speed Cs and 
minimum grid scale. Estimates of the viscosity are trickier, 
but the viscosity for a fully ionized u nmagnetized thermal 
plasma (|Braginskiilli95i : ISpitzei|[l96l) can be written as 

^ rp5/2 

Hum. = 2.2 X 10^^'^ -—--gcm~^ . (3) 



For typical values of the ICM temperature (kT ~10 keV) 
and density and assuming a sound speed of Ci, = 1000 km/s, 
the Reynolds number at our peak grid resolution is roughly 
Re ~ 250. Flows with these Reynolds numbers are below the 
expected threshold for transition from laminar to turbulent 
flow in an ideal gas. Additionally, for our simulations, the 
details of turbulent flow (should it exist) are not captured 
at our smallest grid scale or below. However, Re scales lin- 
early with the characteristic length of the flow, and large 
scale eddies should be captured. We expect that turbulence 
injection in the ICM hap pens at scales well above o ur grid 
size (L ^ 100 - 400fcpc) (|Brunetti fc Lazarian|[201ll '). What 
we do capture in this calculation is the kinetic energy as- 
sociated with bulk random motions of the ICM gas at and 
above the grid scale of the simulation. This motion would 
be converted to thermal energy through the turbulent cas- 
cade (and through shock dissipation) in the real ICM. So 
this calculation is made to simply characterize the available 
energy in random bulk motion of the ICM, which we as- 
sume will eventually be dissipated into heat (and possibly 
CRs) through the turbulent cascade. Indeed, recent work by 
IPaul et al] (|201ll ) indicates that our metric k should be rel- 
atively insensitive to the spatial resolution of the simulation 
grid. 

The quantity k is tracked with time in the history of 
each cluster to determine the change in the TKE. If the 
paradigm that cluster radio halos are powered by turbu- 
lent re-acceleration of CRs is correct, increases in the TKE 
should lead to an increase in the radio brightness. Therefore, 
we can make estimates of the lifetime of radio halos by track- 
ing the TKE evolution. We also study the correlation of the 
X-ray structure measures with the TKE; since the presence 
of a radio halo is observed to be linked to X-ray structure, 
the TKE should likewise be correlated to cluster structure 
if indeed turbulence is connected to the production of radio 
halos. 



3 CONNECTION OF CLUSTER MERGERS, 
STRUCTURE, AND TURBULENCE 

3.1 Distribution of Cluster Structures 

ICassano et al] (l201ol ) show that radio halo clusters can be 
distinguished from non-radio halo clusters based on their 
observed X-ray structure. In particul ar, in term s of th e two 
structure measures we consider, iCassano et al.l (|2010i ) flnd 
that clusters hosting radio halos tend to lie in the upper- 
right quadrant of the w — P3/P0 plane with high centroid 
shifts and high P3/P0. Below, we compare the distribution 
of clusters structu res for our simu l ated clusters to observed 
clusters from the ICassano et al.l (|2010l ) sample and show 
that they are similar. 

In Figure 1, we show the distribution of structures of 
our simulated clusters in the w — P3/P0 plane for all time 
output s with 0.2 ^ z ^ .4, the r edshift range s p anned 
by the [c assano et al.l (|2010D sample. [C assano et al.l (|2010D 
make cuts for radio halo versus non-halo clusters based on 
the median of each structure measure, which for their sample 
gives w = 0.012 and P3/P0 = 1-2 x 10"''. Most of the radio 
halo clusters in their sample do in fact have w > 0.012 and 
P3/P0 > 1-2 X 10"''. For our simulated sample in a similar 
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Figure 1. The centroid shift w versus P3/P0 plane for all the 
simulate d cluster images in th e redshift range 0.2 ^ ^ ^ 0.4 (the 
range m ICassano et alj fe OlO")). Rod lines indicate the structure 
cuts from lCassano^t al.l t2010.) , black are the median values for 
our simulated clusters. 



redshift range (0.2 ^ z 0.4), we find median values of 
w = 0.010 and P3/P0 = 4.6 x 10~^, similar p articularly for 
the centroid shift to the ICassano et al.l (|2010l ) values. Both 
the lCassano et al] (|2010| ) structure cuts (red lines) and the 
median values in the simulations (black lines) are shown in 
Figure 1. We find that in the 0.2 < 2 ^ 0.4 redshift range, 
24 % of the simulated clu ster images have structures above 
the ICassano et al.l (|2010l ') cuts {w > 0.012 and P3/P0 > 
1.2 X 10~^) while 36% have structures above the median 
values for the simulated sample. 

Overall there is a good match between the structure of 
clusters in our simulations and observed clusters in terms 
of both the typical cluster structure and ran ge of structure 
measu res. Comparing Figure 1 to Figure 1 in lCassano et al.l 
l|2010h gives the visual impression that a larger fraction of 
simulated clusters lie in the off-diagonal (lower-right and 
upper-left) regions of the structure plane implying a less 
strong correlation of w and Pa/Po- In Figure 1, a total of 
26 — 28% (depending on the cuts used) of the simulated clus- 
ter images lie in th e two off-diagonal qua drants compared to 
5 out of 32 in the ICassano et al.l (|2010l ) sample, not a sig- 
nificant difference given the sample size. 

The fraction of simulated clusters in the upper-right 
quadrant of the w — P3/P0 plane, which have X-ray struc- 
tures similar to those observed for radio halo clusters, is 
a good match to the 29 ± 9% of X-ray luminous clusters 
found to host radio halos in the survey of a complete cluster 
sample with the Gia nt Meterwave Ra dio Telescope (GMRT ) 
l|Venturi et all 120081 ) (from which the lCassano et all (|2010l ') 
sample is drawn) and t o the 12 out o f 32 c lusters hosting 
giant radio halos in the lCassano et al] (|2010l ) sample. 

If we consider instead the full redshift range of the sim- 
ulations (0.0 ^ 2 ^ 0.9), we find that clusters have ra- 
dio halo like structu res 19% (31%) of the time given the 
ICassano et al.l (|2010l ) (simulated) structure cuts. As we will 
show in later sections, some clusters spend most of their 
time in disturbed states while other clusters are relaxed over 
nearly all of their history. 



3.2 Cluster Histories 

We follow the evolution of both the structure measures and 
the TKE for individual simulated clusters through time. 
This reveals a wide range of cluster histories. We investi- 
gate the effect of mergers on measured cluster structure (i.e. 
the observable effects of mergers) and on the TKE. We can 
then ask when clusters have structure measures similar to 
those observed for radio halos and how this relates to the 
TKE and merger history. 

The time histories of the TKE in 6 of the simulated 
clusters in our sample are shown in Figure [2] The clusters 
are chosen to represent a wide range of time histories. We 
use the kinetic energy density as defined earlier to remove 
the effect of increasing mass on the result. In each case, for 
the purposes of comparison with the structure measures, we 
only examine the value of k inside a radius of 500 kpc. It 
is clear from these histories that the value of k undergoes 
major variations as a function of time in each cluster. These 
variations can change the value by factors of 5 or more over 
time intervals of 0.5-1 Gyr. The strong increases in the TKE 
will be shown in later sections to correspond directly to the 
incidence of merging events. These merging events, as has 
been shown in previous investigations, also correspond to 
changes in the X-ray structure measures. 

However, for the structure measures, the variations with 
time for individual clusters are not as smooth as those for the 
TKE. Figure |3] shows the variation in the P3/P0 power ratio 
as a function of redshift for the same 6 simulated clusters 
shown in Figure (21 While the TKE rises and falls smoothly 
during merger events, the value of the X-ray P3/P0 power 
ratio jumps up and down significantly during the whole life 
of the cluster. From previous studies, we were aware that 
the power ratios can vary significantly on short time scales. 
Merger phase can have a strong effect on the structure mea- 
sures. Strong dips in P3/P0 are often seen, for example, at 
the moment of core passage in major mergers. The X-ray 
structure measures tend to be highest in the pre-merger 
and post-core passage re-expansion phases when subclus- 
ters are well-separated and/or asymmetric features are most 
prominent. Statistically, X-ray structure measures like w 
and P3/P0 correlate strongly with cluster dynamical state, 
but significant scatter is present for individual cluster mea- 
surements (projectio n effects also contribute to the scatter) 
(iJeltema et al. 1 120081 ). Therefore, though a correlation be- 
tween the turbulence kinetic energy and the X-ray structure 
measures is expected, it is not immediately obvious from 
individual cluster time histories that this is the case. 

To further illuminate the effects of cluster mergers, in 
Figures |4l (5] and [6] we show three example mergers from our 
simulated cluster sample which exhibit quite varied merger 
progessions. During the mergers, we have measured the value 
of both (w) and P3/P0, as well as the turbulence kinetic 
energy k. In each figure, the images from left to right in 
both the upper and second row are a time series of X-ray 
surface brightness (upper) and X-ray weighted temperature 
(second). The two plots in the second row from the bottom 
show the progression in time of the TKE (left) and the po- 
sition in the P3/Po-{w) plane. The bottom panels show the 
evolution of the individual structure measures as a function 
of redshift. Letters indicate the position in the time series. 
In Figure |4l we see in panel A a subhalo approaching the 
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main cluster; in panel B, this subcluster is driving an obvious 
large shock into the gas. Point B in the TKE time series is 
associated with the peak of the enhancement, and also the 
structure measures have moved into the disturbed area of 
the P-i/Po-((w) plane. As the subcluster passes through the 
main cluster, the TKE declines gradually, and the struc- 



ture measures move back into the "relaxed" regions of the 
structure plot. Between points B and C, the structure mea- 
sures show the typical progression seen in cluster mergers: 
both P3/P0 and (w) show a peak near point B when the 
subcluster enters the 500 kpc aperture; as the cluster un- 
dergoes core passage between points B and C, the structure 
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Figure 4. Merger snapshots for a simulated galaxy cluster in our sample. Upper panels show images of the synthetic 0.3-8.0 keV X-ray 
surface brightness while the panels immediately below show the projected X-ray weighted temperature. Image panels are roughly 1 kpc 
across. The black circle overlaid on the final time step indicates the fiducial 500 kpc radius used in our analysis for that timestep. The 
panels in the third row show the time evolution of the turbulence kinetic energy (left) and the location of the cluster in the P3/P0 — (ui) 
plane for the four time snapshots as indicated by letter. Bottom pan els show the time evolution of P3/P0 and (ui) with redshift; green 
dashed lines indicate the structure cuts used bv lCassano et all 1I2OIOI) to select radio halo clusters. The mass of the main cluster ranges 
from 4.0x10"Mq at (A) to 5.0x10^^Mq at (D). 



measures both dip; P3/P0 and (w) then show a second peak 
following core passage in the re-expansion phase before de- 
clining again. Point C is near the end of re-expansion and 
outgoing shocks can be seen in the temperature map. Com- 
pared to the structure measures, the TKE evolves smoothly 
throughout the merger. At point D, the small remaining core 
of the subcluster can been seen re-entering the cluster, but 
this more minor event does not lead to as large a rise in the 
structure measures or the TKE. 



Figure \S\ shows a slightly more complicated merger. 
First, we have two approaching subclusters (seen in the up- 
per center and upper right of panel A). The first merger 
drives a large shock, seen in the temperature image of panel 
B. This again leads to a spike in the TKE and a movement 
of the structure measures in to the disturbed regime. Before 
the cluster can relax, the second, apparently smaller subclus- 
ter plunges in, enhancing the TKE again and again boosting 
the structure measures toward the disturbed regime (point 
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Figure 5. Merger snapshots for a simulated galaxy cluster in our sample. Upper panels show images of the synthetic 0.3-8.0 keV X-ray 
surface brightness while the panels immediately below show the projected X-ray weighted temperature. Image panels are roughly 1 kpc 
across. The black circle overlaid on the final time step indicates the fiducial 500 kpc radius used in our analysis for that timestep. The 
panels in the third row show the time evolution of the turbulence kinetic energy (left) and the location of the cluster in the P3/P0 — (w) 
plane for the four time snapshots as indicated by letter. Bottom pan els show the time evolution of P3/P0 and {w) with redshift; green 
dashed lines indicate the structure cuts used bv l Cassano et al. l S to select radio halo clusters. The mass of the main cluster ranges 
from 2.8x10"Mq at (A) to S-SxIOI^Mq at (D). 



C). After an extended period, the TKE declines and the 
structure measures move in to the relaxed quadrant of the 
P3/P0 — (w) plane (point D). 

Figure |S] shows a very complex scenario. In panel A we 
see two subhalos in the center left of the image which are 
in the process of merging. In panel B, the subhalos have 
merged with each other and are in a re-expansion phase 
(note the merger shock in the temperature image of panel 
A, and the outgoing shocks in B). They then pass by the 



main cluster on either side, with the leading shock sweeping 
over the main cluster. Meanwhile, in panel C, note another 
approaching subcluster to the upper right, which by panel 
D has plunged into the main cluster. This multiple merger, 
as seen in the context of the TKE and structure measures, 
has an increase in k during both merger events, while the 
structure measures remain high throughout. This is a very 
dynamic environment, where multiple mergers take place 
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Figure 6. Merger snapshots for a simulated galaxy cluster in our sample. Upper panels show images of the synthetic 0.3-8.0 keV X-ray 
surface brightness while the panels immediately below show the projected X-ray weighted temperature. Image panels are roughly 1 kpc 
across. The black circle overlaid on the final time step indicates the fiducial 500 kpc radius used in our analysis for that timestep. The 
panels in the third row show the time evolution of the turbulence kinetic energy (left) and the location of the cluster in the P3/P0 — (ui) 
plane for the four time snapshots as indicated by letter. Bottom pan els show the time evolution of P3/P0 and (ui) with redshift; green 
dashed lines indicate the structure cuts used bv lCassano et al] l|2010h to select radio halo clusters. The mass of the main cluster ranges 
from a very small 1.6xlOl'*M0 at (A) to more than double its original mass, 3.5xlO^*Af0 at (D). 



throughout the cluster lifetime, leading to long periods of 
elevated TKE and disturbed structure measures. 

In general, we find that there is no "typical" cluster 
history. Some clusters undergo complex mergers or progres- 
sions of successive mergers leading them to be disturbed over 
most of their history, while other clusters are nearly always 
relaxed, experiencing only one or two minor mergers from 
z = 0.9 to the present. This diversity highlights the impor- 



tance of using cosmological scale simulations, since idealized 
mergers would not capture the full complexity of cluster evo- 
lution. 
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3.3 Assessing the Turbulent Re- Acceleration 
Model 

It is clear from the cluster histories that the TKE can in- 
crease drastically during cluster mergers (as you would ex- 
pect) and that clusters can display quite different TKE evo- 
lutionary histories. In this section, we consider how the TKE 
correlates with the observed properties of radio halo clusters 
and if this is consistent with the turbulent re-acceleration 
model of radio halo production. Specifically, we examine the 
TKE in clusters with and without X-ray structures consis- 
tent with those observed for radio halo clusters, and we look 
at the correlation between total TKE and cluster X-ray lu- 
minosity. 



3.3.1 Correlation of Cluster Structure and Turbulence 

One of the main relevant questions of our investigation 
is whether the presence of enhanced TKE correlates with 
the value of the structure measures for the X-ray images. 
We show in Figure [7| the P^/Po-iw) plane for each of the 
clusters shown in Figure [2] Each data point in these im- 
ages represents one time interval from the simulated data 
for the cluster in question. The data points are colored 
by their relative TKE, to indicate the level of turbulence 
at the time the structures are measured. The colors go 
from red {k < 2kmin) to yellow {2kmin < k < Skmin), to 
green (Sfemin < k < 'ikmin), to blue (fc > Akmin). In each 
case, kmin represents the minimum value of k over the full 
time interval for which the cluster and its predecessors have 
M > IO^^Mq. Using kmin may lead to scatter in the results, 
however, it is not clear from any of the histories what the 
base value for k should be. In many objects, increases in k 
due to a new merger occur before k has finished decreas- 
ing from prior mergers. Certainly, by eye the visual effect 
of mergers lasts longer than the time between mergers for 
many of the objects. 

The resulting correlations vary from cluster to cluster, 
but in general, Figure[7]shows an obvious trend for more dis- 
turbed power ratios and centroid shifts to be accompanied 
by a higher TKE state. We can quantify this result further 
by examining the correlation of cluster structure with the 
TKE in a global sense across all the clusters. One key mea- 
surement is the amount of time spent in both an elevated 
TKE state and in the upper right quadrant of the P3/P0 
- (w) plane. Conversely, we can also check that the lower 
left quadrant of this plot is more likely to be populated by 
clusters with lower TKE. 

Figures [8] and [9] show the correlation between the tur- 
bulence kinetic energy and the location of the cluster's pro- 
jected structure measures. These plots show the fraction of 
time for each cluster that it spends in an elevated TKE state 
(either twice or three times the minimum of k) versus the 
time spent in the upper right quadrant or lower left quadrant 
of the P-i/ Po-{w) plane. It is clear that these two measures 
are correlated, though there is significant scatter. To first 
order, the clusters that spend more time in the elevated k 
states spend more time with disturbed structure measures. 
It is also evident that the converse is true, those that spend 
more time in elevated k states spend less time in the lower 
left of the structure plane. It is also obvious from these Fig- 
ures that the range of values for fractions of time spent in 



Table 1. Fraction of Time with Elevated TKE for Location in 
Structure Plane 



Ps/Po-iw) Quad 


2x TKE 


3x TKE 


4x TKE 


Upper Rigiit 


0.92 


0.77 


0.58 


Lower Left 


0.44 


0.19 


0.13 



various regimes are large. Some clusters spend a large frac- 
tion of time with elevated TKE and disturbed morphology, 
while others spend almost no time in these states. It is im- 
portant to note that clusters do not follow some standard 
evolutionary path, where they each spend a typical fraction 
of time in elevated TKE states, but that each evolutionary 
track is unique. There is a wide range of evolutionary histo- 
ries, even in a sample this small. 

Table[T]shows, for the whole sample of clusters and time 
outputs, the fraction of time that the structure measures lie 
in both the relevant quadrants of the Pz/Po-{w) plane and 
in elevated TKE states. The table shows clearly that clusters 
that have disturbed structure measures are much more likely 
to have elevated TKE. While this correlation is not perfect, 
if a cluster has disturbed structure measures, it has a 92% 
chance of having its TKE at least double the minimum over 
the interval in question. Likewise, only a small fraction of 
clusters with structure measures indicating that they are 
relaxed have very elevated TKE. Though relaxed structure 
measures can not perfectly predict low TKE, a number of 
factors contribute to this problem. First, the use of kmin as 
the base level for TKE almost certainly adds scatter to the 
correlation. Also, as we have shown in earlier sections, there 
are moments (particularly at merger core passage) where 
structu re measures look relax ed, yet the TKE is elevated. 

In ICassano et al] l|20ld ). a strong case is made that 
there is a real correlation between incidence of radio ha- 
los and location in the planes described by combinations of 
the concentration parameter, (w) and P3/P0. The correla- 
tion found found here between cluster structure and e l evated 
TKE lends credence to the claims of ICassano et al.1 (l2010l ) 
and others that turbulence could indeed be the energy source 
for the CRs that create the radio halo emission. 



3.3.2 X-ray Luminosity versus Turbulence Kinetic Energy 

Under the assumption that turbulent bulk motions provide 
the energy source for the CR electrons radiating as radio 
halos, we plot the total (not specific) turbulence kinetic en- 
ergy Kt = Mgaak against the synthetic X-ray luminosity of 
our simulated clusters in Figure 1101 The X-ray luminosity 
is extracted from the projected maps used for the struc- 
ture measures, and the turbulence kinetic energy is inte- 
grated in each cluster volume (radius of SOOkpc). The rela- 
tion is shown in Figure[Tn]for a randomly selected 10% of the 
timesteps in each cluster history (to avoid the appearance 
of evolutionary tracks). We can not expect identical behav- 
ior to the Lx — Pi.iGHz relation for observed clusters due 
to the complicating factors of the magnetic field strength 
and distribution, and the aging of the accelerated particles. 
However, it is interesting to check whether the "switch" for 
on/off states of radio halos results directly from a quick tran- 
sition in turbulence state. Previous work (e.g., IPaul et al.l 
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Figure 7. The P-j,/ Pq-{w) plane for projections of all time intervals where M200 > 10^* Mq for the six 2=0 clusters whose TKE time 
histories a re shown in Figured In this plot, [w) is the centroid shift with the core excluded. Lines indicate the fiducial cut lines in P3 /Pq 
and [w) in lCassano et ahl 1I2OIOI) . Points are colored to indicate the relative level of turbulence kinetic energy for each time interval. Red 
indicates intervals where k < 2kmin, where fcmin is the minimum value of the TKE over the full time range for a given cluster. Yellow 
indicates 2kmin < k < Sfcmim green is Sfcmin < k < Ak^in, and blue is for k > Akmi„. 
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Figure 8. Correlation of the time spent by our simulated clusters 
in elevated turbulence kinetic energy states with the fraction of 
time spent in the upper right (disturbed) qua drant of the P3/P0- 
(w) plane (as defined bv lCassano et aU 1I2OIOI) . Green points rep- 
resent a value of A; > 2kmin: blue indicates k > Skmin- Note that 
some clusters spend almost the entire interval we studied in a 
high TKE state, while others spend almost no time at all in that 
state. 

I2OIII ) suggests that turbulent pressure is long-lived in the 
ICM after mergers. Our relation shows no obvious evidence 
for bimodality in the relation between X-ray luminosity and 
turbulence kinetic energy. Indeed, our result suggests that 
the transition from turbulent to relaxed is gradual, as in- 



Figure 9. Correlation of the time spent by our simulated clusters 
in elevated turbulence kinetic energy states with the amount of 
time spent in the lower left (relaxed) quadrant of the P:j/Po-{w) 
plane (as defined bv .Cassano et al.it201Q) . Green points represent 
a value of fc > 2kmin, blue indicates k > 'ik^in- 



dicated by our earlier time histories. Our analysis shows 
that there is an order of magnitude spread in the value of 
Kt f or a given X-ray lumino s ity. T his result is consistent 
with iRudnick fc LemmermanI l|2009l ) , which indicates that 
the value of Pi. 4 can take on a range of values for a given 
Lx, under some upper envelope. It is possible that once a 
large population of radio halos can be detected at longer 
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Figure 10. Synthetic 0.3-8.0 keV X-ray luminosity versus tiie 
total turbulence kinetic energy (Mgask) for a randomly selected 
10% of the time outputs for all clusters and their progenitors in 
our sample. Lx is determined from the projected mock images, 
Kt is extracted from the three-dimensional simulation data. 



Figure 11. Histogram of all continuous time intervals for clusters 
in our sample where the turbulence kinetic energy inside r=500 
kpc is elevated to A; > 2kmin- The median interval length is 1.12 
Gyr. 



wavelengths (given their steep radio spectra) any clear bi- 
modality in the L^-Pr relations will be washed out. Also, 
though the Lx — Kt relation presented here is not bimodal, 
there are other relevant physical efTects on the radio emis- 
sion (magnetic field strength and variation, presence and 
spectrum of a pre-existing electron population) which are 
not included here that could induce strong effects on the 
resulting Lx — Pr ratio. For example, Brunctti ct al] (|2009l ) 
suggest that as the turbulent kinetic energy decays following 
a merger, the the suppression of the synchrotron emissivity 
in a fixed radio band will be strong. This also suggests that a 
factor of two difference in the value of k will lead to a factor 
of 10 suppression in the synchrotron luminosity at 1.4 GHz. 
In that case, a continuous distribution in Kt for a given Lx, 
as appears in our figure, will still show a sharp boundary, 
leading to the observed bimodality. 



4 DISTRIBUTION AND DURATION OF 

TURBULENCE IN SIMULATED CLUSTERS 

If in fact turbulent re-acceleration is responsible for generat- 
ing the relativistic electrons powering radio halos, then the 
evolution of the TKE in our simulated clusters can tell us 
about the possible lifetime of radio halos and the frequency 
with which clusters are expected to host radio halos. In this 
section, we examine the typical duration of episodes of ele- 
vated turbulence (i.e. the merger relaxation time) and the 
fraction of clusters at a given time that are in elevated TKE 
states. 

4.1 Lifetime of Merging Events 

We have run an analysis to determine the typical timescale 
for clusters in our sample to be in elevated TKE states. We 
take every cluster's time history and identify all time in- 
tervals over which the TKE is continuously elevated to at 
least 2kmin. The histogram of those intervals is shown in 



Figure [TT] The median length of time that a cluster spends 
in a continuously elevated TKE state by this definition is 
tele ~ 1.13 Gyr. If we redefine the intervals by elevation 
to Skmin or above, then teie = 0.77 Gyr. For k > Akmin, 
tele = 0.77 Gyr. While it may not be obvious why the me- 
dian interval for the highest two k bins should be the same, 
the rough shape of the distribution is preserved at these 
ranges, though the normalization has declined. So it is safe 
to say from our sample that the median time interval over 
which the TKE is elevated inside r =500 kpc is around 1 
Gyr, or slightly less depending on the threshold used. How- 
ever, Figure [TT] also shows a tail to very long intervals of 
elevated TKE with some clusters having time periods of 5 
or more Gyrs with continuously elevated TKE. This result 
again shows the diversity of cluster histories: in the turbu- 
lent re-acceleration model some clusters would be expected 
to host radio halos over most of their lifetime while other 
clusters never would. 

Continuous time intervals for elevated structure mea- 
sures show a different result. T he median value for inter- 
vals of P3/P0 elevated above the lCassano et all (|2010t ) limit 
is much shorter than the elevated k intervals at teie ~ 
O.SSGyr. For the centroid shift we see a similar result of 
tele ~ 0.39Gyr. This result is not particularly surprising 
given the time histories shown in Figures [l] [5] and (6) The 
values of the structure measures are more highly variable 
over short time scales than the value of k, and can show 
strong variations at core passage for instance. If we identify 
by eye time periods with low structure measures associated 
to core passage events and interpolate over them, a longer 
typical timescale for elevated structure measure of around 
O.bGyr is found, though this procedure is subjective. 

4.2 Fraction of Clusters with Elevated Turbulence 
Kinetic Energy 

If the energy source for giant radio halos is turbulent re- 
acceleration, an interesting measure is the fraction of ob- 
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Figure 12. The fraction of simulated massive (M > 5 X 10^* Mq) 
clusters as f{z) that have k > 3kmin- Shaded region shows Pois- 
son error bars due to the finite sample of simulated clusters at 
each redshift. 

jects with elevated turbulence kinetic energy. In Figure [121 
we show the fraction of massive clusters (M^ 5 x 10^"^ Mq) 
as a function of redshift that have k elevated by an integer 
factor above the minimum over their history. The fraction is 
a noisy, but gradually declining function of time depending 
on what threshold of k is used. For k > 3kmin as shown 
in Figure 1121 the fraction of clusters with elevated TKE at 
high redshift {z > 0.4) is noisy, but high (50-100%), and 
declines to roughly 20% at z = 0. The trend for values of 
k > 2kmin and k > Akmin is similar, but the normaliza- 
tion changes. For k > 2kmin the typical fraction of clusters 
declines from 60-100% at the higher redshifts to 35-50% at 
z=0. For k > 4:kmin the fraction is 50-100% at high z, with 
a drop to around 8% at z=0. Our results imply that the 
fraction of massive clusters that are energetically favorable 
for producing CRs through turbulent re-acceleration is de- 
clining with redshift. 

Observational searches find a fraction of roughly 30% of 
cluste rs at moderate reds hifts (0.2 < z < 0.4) hosting radio 
halos (jVenturi et al ]|2008h . This result will obviously depend 
on the survey depth and should also depend strongly on the 
frequency of observation, as steep spectrum radio sources are 
much brighter at longer wavelengths. In addition, we should 
point out that radiative losses to relativistic particles due to 
inverse Compton (IC) interactions with the CMB scale as 
(1-1- 2;)*, strongly limiting the possibility of radio hal os at the 
highe r redshifts (Cassano & Brunctti 2005; Cass ano et al.l 
I2OO6I ). Given the uncertainty in the energetic threshold for 
TKE to produce a radio halo, our range for the fraction of 
clusters with elevated TKE in the redshift range of most 
observed radio halos (0.2 < z < 0.4) of 20-50% is broadly 
consistent with the observed ratio. 

5 DISCUSSION 

Observations show a connection between the detection of 
diffuse, large-scale radio halo e mission and t he presence of 
merging activity in clusters (e.g. lBuot3l200ll : ICassano et al.l 



[20TO). In particular, these studies show a correlation be- 
tween observed X-ray cluster structure, quantified based on 
the centroid shift and the power ratios, and the presence 
of radio halo emission. While the origin of radio halos is 
still a matter of debate, it has been suggested that merger 
generated turbulence could accelerate a pre-existing rela- 
tivistic electron population (from supernovae, AGN, previ- 
ous merger/shock activity) to energies sufficient to produce 
the radio emission. 

In this paper, we study the evolution with time of mas- 
sive clusters formed within a numerical hydrodynamic, cos- 
mological simulation and investigate whether a subset of 
the X-ray and radio observational evidence from clusters is 
consistent with the turbulent re-acceleration model for ra- 
dio halos. In particular, we track the evolution of random 
bulk ICM motions, the kinetic energy of which is dissipated 
through shocks and the turbulent cascade. We study the 
relationship between turbulence kinetic energy and X-ray 
structure measures, how these evolve during mergers, and 
whether the level of turbulence in simulated clusters with 
radio halo-like cluster structures supports the turbulent re- 
acceleration model. 

We show that the distribution of the X-ray structure 
measures (cent roid shift w and th ird order power ratio 
P3/P0) used bv lCassano et al.l (|2010l ') to separate radio halo 
and non-radio halo clusters are reproduced fairly well in our 
simulated cluster sample and that a similar fraction of sim- 
ulated clusters have "disturbed" X-ray structures compared 
to observations. We also find that clusters in the disturbed 
regimes of the structure plane are highly likely to have el- 
evated turbulence kinetic energy (at the scales we resolve). 
For some clusters, we have shown this graphically (Figure 
(Tjl. Quantitatively, clusters with structure measures in the 
"disturbed" part of the structure plane have a 92% chance 
of having an elevated k at least two times above the mini- 
mum over the cluster lifetime. Though individually, there is 
a significant amount of variation in the structure measures 
for a given va l ue of fc, as a group they correlate strongly. 
ICassano et alj l|2010l ) suggest that the fact that radio halo 
clusters are more disturbed in terms of their structure mea- 
sures results from a difference in the amount of turbulence 
in disturbed versus relaxed clusters and that turbulent re- 
acceleration of CR electrons is the source of radio halos 
in galaxy clusters. The strong correlation between cluster 
structure and elevated turbulence in simulated clusters sup- 
ports this theory. 

The value of the turbulence kinetic energy k undergoes 
large variations during mergers of factors of a few to several. 
Mergers are typically well-defined in k{z), with k smoothly 
increasing and then decreasing through the merger, though 
additional mergers can (not infrequently) occur before the 
system can fully relax. The time evolution of the structure 
measures during mergers are not as well defined and are 
a strong function of merger phase. Both P3/P0 and (w) axe 
noisier than fc as a function of time in our simulated clusters. 
A classic example in many of our simulated clusters is the 
sharp reduction in the values of the structure measures as 
a merger proceeds through the core passage of a subcluster. 
Though the turbulence kinetic energy is high at this stage, 
the object could be called "relaxed" using the structure mea- 
sures and significant substructure is not visible in the X-ray 
image. The result of this difference is that the time scales for 
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continuously elevated k are of order 2-3 times longer than 
those for elevated structure measures. The median interval 
for enhanced k to twice the minimum for a given cluster is 
approximately t—1 Gyr. For Pg/Po and (w), this number is 
closer to r= 0.33 Gyr. For these reasons, translating from 
the observed structure measures to the dynamical state, or 
more specifically to the energetics of the bulk motion of the 
gas, has significant scatter for individual clusters. However, 
statistically these properties correlate well. 

Though prior studies suggest tu rbulent kinetic en ergy 
is very long-lived in clusters (e.g., IPaul et al] 1201 ll ). we 
show that k drops significantly in cluster volumes over much 
shorter timescales (~ 0.5 Gyr ). However , these results may 
be consistent, given that the IPaul et al] (|201lh result uses 
a different criterion (turbulent pressure > 5% of the total 
pressure) than we do for elevated turbulence. The relatively 
quick drop in TKE, coupled to short synchrotron cooling 
times, may account for the sharp difference in X-ray struc- 
ture for RH versus non-RH clusters. 

An interesting result of this work is that we find that 
many clusters spend a large fraction (~80-|-%) of their 
time in disturbed states, while others spend very little time 
(~10%) in such states. Thus, it may be the case that some 
galaxy clusters are rarely, if ever radio halo clusters during 
their lifetime, while others almost always host radio halos. 
While the typical timescale for the turbulence kinetic en- 
ergy to be elevated is 1 Gyr, these periods can be as long 
as 5 or more Gyrs in cluster which experience complex or 
repeated mergers. 

Finally, we find no obvious evidence for bimodality of 
the Lx-K relation (where K could be a proxy for radio 
power) in our simulated sample of clusters. However, this 
is not necessarily surprising given that there are a number 
of reasons why such a gap may exist in the observed L^-Pia 
relation that are unrelated to the dynamical state of the 
gas. For example, we have not modeled magnetic field varia- 
tions, and the turbulent re-acceleration mechanism requires 
a pre-existing population of CRs on which t o operate, both 
of wh ich may vary from cluster to cluster. iBrunetti et al.l 
l|2009t ) have shown that the relationship between the kinetic 
energy in turbulence and the fixed band radio luminosity 
is not linear. Their work suggests that a relatively small 
drop in TKE results in a large change in the synchrotron 
luminosity, which could easily account for the observed bi- 
modality, yet still be consistent with our result here. To fully 
understand this issue, it is necessary to follow the magnetic 
fields and CR particles explicitly in the simulation to see 
the true evolution of the radio emission. A full treatment 
of simulating a large sample of galaxy clusters and their ra- 
dio emission must handle the injection of CRs into the ICM 
from reasonable sources like AGN, SNe, and large-scale ac- 
cretion shocks. This work we defer to upcoming simulations 
using a new MHD-I-CR code. 
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